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Abstract. We present adaptive mesh refinement (AMR) hydrodynamical simulations of the interaction between Type la 
supernovae and their companion stars within the context of the single-degenerate model. Results for 3D red-giant companions 
without binary evolution agree with previous 2D results by Marietta et al. We also consider evolved helium-star companions 
in 2D. For a range of helium-star masses and initial binary separations, we examine the mass unbound by the interaction and 
the kick velocity delivered to the companion star. We find that unbound mass versus separation obeys a power law with index 
between -3.1 and -4.0, consistent with previous results for hydrogen-rich companions. Kick velocity also obeys a power-law 
relationship with binary separation, but the slope differs from those found for hydrogen-rich companions. Assuming accretion 
via Roche-lobe overflow, we find that the unbound helium mass is consistent with observational limits. Ablation (shock 
heating) appears to be more important in removing gas from helium-star companions than from hydrogen-rich ones, though 
stripping (momentum transfer) dominates in both cases. 
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INTRODUCTION 

Type la supernovae (SNe la), unlike core-collapse supernovae, are most likely purely a phenomenon of close binary 
systems. Their defining absence of hydrogen spectral lines and presence of silicon lines suggest an origin in compact 
stellar remnants, while the progenitor delay times inferred from their presence in early- and late-type galaxies and 
their absence from interarm regions require these remnants to be low-mass, ie. white dwarfs (WDs). The relative 
homogeneity of the SN la population suggests that most white dwarfs explode at a single mass, probably the 
Chandrasekhar mass, 1.44 M [1]. Because the distribution of white dwarfs peaks at 0.6 M Q [2], Chandrasekhar- 
mass models require that the white dwarf accrete from a binary companion. 

Most work on progenitor models for SNe la has therefore focused on two general cases: the single-degenerate model, 
in which a white dwarf accretes matter from a normal stellar binary companion, then becomes unstable to explosive 
nuclear burning [1, 3], and the double-degenerate model, in which two white dwarfs merge due to gravitational- 
wave emission [4, 5]. Single-degenerate models must cope with the problem of the companion star's hydrogen and 
helium; observational limits on H/He abundances in SNe la are very constraining [6, 7]. These models also have the 
difficulty that the mass accretion rate onto the white dwarf must be in a fairly narrow range above 10~ 7 M yr~' to 
avoid either nova outbursts or development of a white-dwarf wind that limits the accretion efficiency [3, 8, 9]. This 
makes it difficult to explain the observed low-redshift SN la rate of 3 x 1CT 5 Mpc~ 3 yr~' [10] using single-degenerate 
progenitors alone. The predicted and observed abundances of double-degenerate systems are sufficient to explain the 
observed rate [11, 12, 13]. However, such collisions may lead to production of an O-Ne-Mg white dwarf followed by 
accretion-induced collapse to a neutron star [14, 9, 15, 16]. It is possible that multiple progenitor channels contribute. 
The observed redshift distribution of SNe la is consistent with at least two populations having different delay times 
measured from the zero-age main sequence [10]. 

We are examining the single-degenerate model to determine, among other things, whether the problem of H/He in 
the companion envelope can be overcome. Single-degenerate progenitors can explode at the Chandrasekhar mass for 
main-sequence (MS), red-giant (RG), or helium-star companions. If the companion is a helium star, the WD may also 
rarely explode below the Chandrasekhar mass via explosive burning of an accreted helium layer that ignites interior 
carbon burning [17]. The Chandrasekhar-mass helium-star channel is regarded as a candidate for the short delay-time 
population of SNe la because this case requires a massive initial MS star [18, 19, 20, 21], while the other two types of 



nondegenerate companion result from systems with low-mass stars and thus should produce longer delay times. 

Multidimensional hydrodynamics simulations of the different channels are difficult because of the large dynamic 
ranges required and the large binary parameter space to explore. Simulations following the ejecta evolution within 
the progenitor system starting from a pointlike explosion potentially make the most direct contact with supernova 
observations and progenitor searches, though they do not address the question of how the supernova explosion begins. 
Such simulations can determine the amount of mass that should be unbound from the companion, determine whether 
hydrogen lines should be seen, determine the level of asymmetry introduced in the supernova remnant by interaction 
with the companion, and provide clues regarding the appearance of the post-supernova companion. 

Marietta et al. [22] explored the RG, subgiant, and MS single-degenerate channels using 2D Eulerian hydrodynamics 
simulations. They found that significant quantities of hydrogen should be unbound from the companion star envelope 
in each case (15% for MS and subgiant cases, and 98% for the RG), in conflict with observational upper limits [6, 7]. 
They also found that the companion star blocks about 7 - 12% of the SN ejecta, introducing anisotropy in the supernova 
remnant and possibly having consequences for Si II spectral line shapes, but that the kick delivered to the companion 
is smaller than the orbital speed in each case and negligible for the RG core. Fallback of part of the MS or subgiant 
companion's envelope produces a dramatic increase in the companion's luminosity, but relatively little contamination 
by SN ejecta occurs. The remnant in the RG case is a helium pre-white dwarf with an extended hydrogen-rich envelope 
of very low mass (about 0.02 M©) that should appear as an underluminous O/B star for up to 10 6 yr. 

More recently, Pakmor et al. [23] updated these results to consider the effect that pre-supernova binary evolution 
has on the structure of MS companions, as suggested by [24], using 3D smoothed particle hydrodynamics (SPH) 
simulations. Pakmor et al. used initial conditions derived from stellar models [9] that included binary mass transfer 
and thus yielded a more compact MS companion. Because of this, they found a tenfold reduction in the amount 
of unbound mass, bringing the predicted amount back into agreement with observational upper limits. They also 
found strong dependences of the unbound mass M un b OU nd an d kick velocity Vkick on initial binary separation a, with 
^unbound 06 a~ 3 ' 49 and vack a~ 1,45 . Their use of SPH with a relatively small number of particles (~ 10 6 ) caused 
them to fail to reproduce the detailed fluid instabilities seen by Marietta et al., but they performed a convergence study 
using the same MS model as Marietta et al. and found good agreement for the unbound mass and kick velocity. These 
quantities therefore should not be too sensitive to the method used, though we note that the degree of contamination 
of the companion by SN ejecta should depend significantly on the treatment of small-scale instabilities. 

We are revisiting this problem using modern adaptive mesh refinement (AMR) techniques, larger computers, and 
companion models that incorporate the effects of binary evolution. These advances permit us to consider larger 3D 
spatial dynamic ranges using the same type of Eulerian shock-capturing hydrodynamics methods as Marietta et al. We 
also examine the helium-star channel, which has not heretofore been considered using hydrodynamical simulations. 
We have previously reported on 2D simulations of the helium-star channel in [25]; here we summarize these results 
and present preliminary 3D results for the RG channel. 

NUMERICAL METHODS 

We use two separate codes to create progenitor stellar models and to evolve SN explosions within binary systems. 
To construct ID stellar models including mass loss due to binary interactions, we use EZ [26], which is based on the 
STARS stellar evolution code originally developed by Eggleton [27, 28]. We have modified the code to include a mass 
loss term to account for binary mass transfer and wind losses. Using parameter values chosen from [9] (for MS and 
RG companions) or [21] (for helium-star companions), we evolve each companion star model up to the point at which 
the supernova explosion is expected to take place. We then interpolate the ID profiles of density, temperature, and 
isotopic abundances onto the 2D or 3D grid used in the next stage of the calculations. We allow the multidimensional 
stellar model to relax for several dynamical times in order to reduce interpolation errors. 

To simulate the explosion, we employ FLASH 3.2 [29, 30], a parallel Eulerian AMR hydrodynamical code based 
on the piecewise parabolic method or PPM [31] for hydrodynamics and a direct multigrid algorithm [32] for gravity. 
We use a nonideal equation of state (EOS) appropriate to stellar material that tabulates the Helmholtz free energy as 
a function of density, temperature, and composition [33]. We separately track hydrogen, helium, oxygen, and carbon 
outside the SN ejecta, while the ejecta themselves are assigned a pure nickel abundance. 

FLASH uses the block-structured AMR package PARAMESH [34], which manages an oct-tree mesh that is 
distributed in parallel using a space-filling curve. To capture flow discontinuities, we refine blocks based on the second 
derivative of the gas density and pressure. We ignore the second-derivative criterion in blocks where the maximum 
density falls below a preset threshold (10 -6 g cm -3 ) in order to avoid excessively refining low-density regions. To 



cope with features peculiar to each of our models, we employ additional refinement criteria as described below. 

Even with AMR, degenerate RG cores are too small for us to resolve directly with our mesh while still achieving 
high resolution elsewhere in our computational domain. Therefore we replace the core gas in RG companions with a 
spherical cloud of about 2 x 10 s particles with radius three times the smallest zone spacing. The gravitational force on 
and due to the particles is determined using the particle-mesh method [35], with an important modification: all of the 
particles in the cloud move rigidly together with the cloud's center of mass. This arrangement ensures that the mapping 
of particle densities onto the mesh and of gravitational forces onto the cloud's center of mass are free of Cartesian grid 
effects introduced by the use of a cloud-in-cell (CIC) mapping kernel. In hydrostatic equilibrium tests of a single RG 
star, this technique allows us to maintain gas velocities in the envelope less than about 1% of the sound speed for more 
than a dynamical time. To ensure maximum mesh resolution of the core, we force refinement of all mesh blocks that 
contain any particles in the cloud. These methods have proven successful in related simulations of common-envelope 
evolution [36]. The MS and helium star companions do not have degenerate cores, so we do not use particle clouds 
with them, but since their outer density profiles are much steeper, we force refinement near their surfaces. 

After allowing the companion star model to relax on a multidimensional AMR mesh, we introduce an SN la by 
adding a small spherical cloud of high-energy gas. To reduce Cartesian mesh effects, we force the mesh to refine about 
the explosion by two extra levels beyond the rest of the computational volume within a spherical region of radius five 
times the initial radius of the explosion. At this level of refinement the initial explosion radius is about three zones 
across. For Chandrasekhar-mass explosions, we use the W7 model [37] to set the ejecta mass, kinetic energy, and 
thermal energy. The W7 model involves a central carbon deflagration that provides a good fit to observed SN la light 
curves and spectra. In this model the explosion energy is 1.233 x 10 51 erg, and the average radially-directed ejecta 
velocity is 8.527 x 10 3 km s _1 . For our purposes, we assume the entire ejecta mass (1.378 M Q ) consists of 56 Ni. 

RESULTS 

We have created a preliminary 3D simulation of the RG case (Figure 1) for comparison with [22]. The RG in this 
simulation had a 0.36 M core and a total mass of 1.05 M Q . The initial binary separation was 6.2 x 10 12 cm, and the 
minimum zone spacing was about 7.3 x 10 10 cm within a computational volume of size 5 AU (7.5 x 10 13 cm). These 
initial conditions resulted in Roche-lobe overflow, so the system was followed for one orbit with the unexploded WD 
represented by a particle cloud (as with the RG core) to develop a realistic aspherical shape. To initiate the explosion, 
we removed the WD particle cloud and replaced it with a gas sphere as described in the previous section. 

Figure 1 shows the double shock structure that forms in this simulation when the SN ejecta begin to interact with the 
red giant. The forward and reverse shocks are separated by a contact discontinuity. When the reverse shock reaches the 
companion, a bow shock structure forms. We reproduce many features of the RG companion case in [22], but in our 
simulation an emptier central region is left behind by the ejecta, because Marietta et al. imposed their explosion at the 
grid boundary, whereas ours is imposed inside the simulation box. Moreover, because we include the orbital motion 
in 3D, the ejecta become somewhat asymmetric in our calculation. We use the kinetic, thermal, and potential energy 
of the gas in each zone to determine the total amount of bound and unbound gas. In agreement with [22], we find that 
almost all the RG envelope is unbound. (Note that neither simulation included the effect of binary mass transfer on the 
initial companion model.) However, a small amount of hydrogen is still left around the RG core. 

Our investigation of the helium-star channel [25] has been more thorough than for the RG case, but our simulations 
to date have been 2D and thus have not included orbital motion. We assumed axisymmetry and placed the supernova 
and the helium star on the symmetry axis. Because helium stars are much more compact than red giants, for these 
simulations we used box sizes of 1 — 5 x 10 11 cm and minimum zone spacings of 5 x 10 7 cm to 1 .2 x 10 8 cm. 

Qualitatively, our results are similar to the results of [22] and [23] for the MS case, but with a more compact com- 
panion and smaller binary separation. We performed a parameter survey using four different helium-star companions 
suggested by [21], described in Table 1, and varying initial binary separations. The helium stars all have relatively low 
masses at the time of the supernova explosion, but they evolve from MS stars of masses 5 — 8 M & . The white dwarfs 
in these models are the remnants of MS stars with initial masses 2 — 6.5 M Q and have accreted material up to the 
Chandrasekhar mass in the course of the pre-supernova binary evolution. 

Figure 2 shows the unbound mass (defined as total unbound helium) and the helium-star kick velocity (defined as 
the center-of-mass velocity of the bound helium) as functions of binary separation in our simulations. We include 14% 
error bars based on results of a convergence test. The unbound mass can be fit by the relation 



M unbound = C ub a m " b M Q , 




FIGURE 1. Left panel: Isodensity surfaces (left) and particle positions (right) for an RG-WD binary after 3.85 x 10 6 s of 
binary evolution (top) and 12,670 s after a supernova explosion (bottom). Passive tracer particles are blue, while active particles 
representing the WD companion are red. The box size is about 2 AU (3 x 10 13 cm). Right panel: Density distribution in the orbital 
plane for RG model at different times after a supernova explosion. The circle shows the initial RG size and position. Each frame 
spans a distance of 1 AU (1.5 x 10 13 cm). The color bar shows log density in g cm~ 3 . 



where a is the orbital separation, m u \, is the power-law index, and the constant C u b depends only on the helium- 
star model (see Table 1). For comparison, we also plot the power-law relation with index —3.49 found by [23] and 
the data from [22] for MS companions (consistent with an index of —3.14). The power-law indices for our helium- 
star companions vary in a small range and bracket their results, suggesting that the index may be insensitive to the 
evolutionary state of the companion. The normalization of the above relation does appear to be sensitive to the nature 
of the companion star. 

At early times, the helium star velocity varies dramatically, but after about 1000 seconds it settles down to a roughly 
constant value. We use an appropriate time average of this value as the kick velocity. For initial binary separations 
larger than 4 times the helium-star radius /?He the kick velocity could not be adequately determined because of the 
initial potential perturbation of the WD. As obtained by [23] and [24], a power-law relation is also found in our 
simulation and can be fitted by the relation 

Vkick = C kick a mkK *, 

where vyck is the kick velocity, m^ck is the power-law index, and the constant Ckick depends only on the helium-star 
model (see Table 1). However, unlike the situation for the final unbound mass, the slope is very different from that 



TABLE 1. Helium-star models considered in [25] 
with best-fit power-law slopes. 
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FIGURE 2. Unbound mass and kick velocity versus initial binary separation for our different 2D helium-star models. Separation 
is scaled to the helium-star radius Ru e . 



found in the case of MS companions (m^ck = —1-45 in [23], and my^ = —1.26 in [22]). Based on a comparison 
with the analytical model of [24], we argue that this difference in slopes may arise because ablation (shock heating) 
and stripping (momentum transfer) have differing levels of importance in the helium-star case versus the RG and MS 
cases. The relative contributions of these processes to mass loss from the companion can be determined by examining 
the kinetic-to-thermal energy ratio of unbound helium as in [25] or, more accurately, by using passive tracer particles 
(Figure 1) and determining this energy ratio at the instant they become unbound. 

The kick velocities for helium-star companions differ by a factor of 2 or more from RG and MS companions 
when considering separations characteristic of Roche-lobe overflow (arlof ~ 3/?He for the models considered here). 
However, for helium stars the final velocity of the companion is dominated by the pre-SN orbital velocity. Assuming 
the orbital velocity to be perpendicular to the kick, we find that the kick velocity is only about 20% of the resultant. 



DISCUSSION AND CONCLUSIONS 

While hydrogen-rich and helium-rich companion stars show many of the same qualitative effects when struck by 
supernova ejecta, the degree of compactness of the companion plays a very important role in determining how much 
mass is lost from the companion and how much of a kick it receives. It remains to be seen whether the ratio of ablated 
mass to stripped mass also depends on the degree of compactness, as this ratio has not yet been computed for the RG 
and MS cases. The companion's compactness increases when we consider stars that have already lost mass due to 
binary interactions. This suggests that an understanding of the outcomes of common-envelope evolution is crucial for 
identifying the types of companions most likely to encounter a single-degenerate SN la situation. 

Our next steps are to revisit the RG and MS channels via 3D AMR simulations that include the effects of binary 
evolution and to examine the helium-star channel in 3D. 3D simulations give us better estimates of the amount of SN 
ejecta material that remains bound to the companion star, since the development of the instabilities that drive mixing 
depends on grid dimensionality. They should also permit us to determine more accurately the fraction of sightlines to 
the SN that can be obscured by the companion. Using tracer particles in 3D, we can also more accurately determine the 
relative amounts of stripped and ablated material (discussed using tracer fluids in [25]). Finally, by including radiation 
transport we expect to determine what fraction of the hydrogen and helium unbound by the supernova is ionized by 
radiation from the companion star. 
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